Calculation of the number of partitions with constraints on the fragment size 
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This article introduces recursive relations allowing the 
calculation of the number of partitions with constraints on 
the minimum and/or on the maximum fragment size. 
PACS numbers: 02.10.De, 05.10.-a, 64.60.Ak 



I. INTRODUCTION 

A partition is an ensemble of positive integers with 
a given sumQ. A partition can also be seen as a way 
to break a piece of discret matter into fragments. The 
number of partitions of a given integer is a quantity 
which is useful in various fields. In the so-called min- 
imal information model jjl for example, all the par- 
titions have the same probability^] P(n) = l/N(S). 
This result is obtained by application of the minimal 
information principle (or maximum entropy), infor- 
mation being defined as ^ n P(n) lnP(n). When this 
equation is differentiated under the only constraint 
s n s = S, all probabilities are found to be equal 



A. Physical Interest 

In many cases, it is interesting to distinguish one or 
more classes of fragments according to their size. In 
percolation, for instance, the sub-critical events are 
defined by the fact that they contain one particular 
fragment which is refered to as infinite or percolative 
in the sense that it connects the two extremes of the 
lattice ||. In the case of conductive bonds, it allows 
an electric current to circulate between two electrodes 
placed on opposite surfaces of the lattice. In the case 
of a coffee-machine, the percolative cluster defines a 
path which allows the vapor to traverse the grounds. 
In Nuclear Physics, other classes of fragments are dis- 
tinguished. The intermediate mass fragments (imf) 
are ions resulting from the violent fragmentation of a 
composite nucleus (produced by the collision of two 



For example, the integer 3 has 3 partitions : {3}, {2,1} 
and {1,1,1}. 

2 In the following, the total number of partitions of the 
integer S will be noted N(S), a partition will be noted 
n: (m, . . . , ns), n 3 being the number of integers s in the 
partition. 



atomic nuclei). When the nucleus is weakly excited, 
it "evaporates" some light particles so that, at the end 
of the process, we are left with light fragments and a 
heavy evaporation residue which contains almost all 
the charge of the initial nucleus. At very high excita- 
tion energies, the nucleus is completely "vaporized" 
into light fragments. In the energy range between 
these two extremes, the nucleus undergoes "multifrag- 
mcntation" into a large number of fragments of all 
sizes. The multifragmentation process is thus char- 
acterized by the production of IMF. Some multifrag- 
mentation models consider the light fragments (with 
charge less than or equal to 2) as nuclear matter in the 
gazeous phase whereas the evaporation residue or the 
IMF form the liquid phase Q. This terminology is also 
used in the field of phase transition models. Thus, in 
various fields, classes of fragments are defined by their 
size. Hence it is of interest to enumerate the number 
of partitions wit h constraints on the fragment size. 
In Section I B some techniques used to calculate the 



number of partitions of an integer without conditions 
on the size of the fragments will be reviewed. It will 
then be shown how the partitions with const raint s on 
the maximum size of the fragments (Section II A), on 
the minimum size of the frag ments (Section II B) and 
on a size range (Section [jXCj) can be enumerated. 



B. The number of unconstrained partitions 

The total number of partitions of the integer S is 
given approximatively by the Ramanujan-Hardy Q 
formula which leading term is : 
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As can be seen, the number of partitions increases 
very rapidly with S. The exact value of the number of 
partitions can be obtained using one of the following 
recursive formulae (the number of partitions of S into 
M fragments is noted N(S, M) and M is referred to 
as the multiplicity) : 

N(S, M) = N{S - 1, M - 1) + N(S - M, M) (2) 
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m— 1 



From this relation || , the equation giving the total 
number of partitions can be deduced : 
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S Int(S/M)-l 

N(S) = 1+J2 N(S-kM-l,M-l) 



M=2 k=0 



(4) 



The Euler M recursive relation leads to the same 
result : 

N(s) = j2(-v k+1 (ms- s A r J1 ) 

k=l 
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II. NUMBER OF PARTITIONS WITH 
CONSTRAINTS ON THE SIZE OF THE 
FRAGMENTS 

A. Constraint on the maximum size 

The number of partitions of an integer S into M 
fragments with size less than or equal to s max will be 
noted Smax 7V(S', M). It is obtained using a modified 
version of the recursive relation (|^) : 

fl °» N(S, M) N(S — 1,M— 1) 



<N(S, M) 



™*N{S - M 
if S< 
« N(M (s ma: 
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5, M). 



The boundary condition is : 

«m«JV(0,l) = 1. 



(6) 
(7) 

(8) 



These relations lead to the calculation of Smax N for 
any value of S and M. The three right hand terms 
in Eq. (Q) are explained as follows. The partition 
ensemble can be shared into two sub-groups. The first 
one contains all the partitions including at least one 
size 1 fragment. One of these size 1 fragments can 
be removed from each partition in the sub-group. It 
follows that the number of partitions in the sub-group 
can be written as Sia ^ x N(S — 1,M— 1). The second 
group includes the partitions with no size 1 fragment. 
Hence one unit can be removed from each fragment 
without modifying the multiplicity. In the absence 
of any condition on the maximum size, the number of 
partitions in the second group would be N(S—M, M). 
However, among the partitions into M fragments of 
the integer S — M, some have one or more fragments 
with size s max . It is not possible to add one unit 
to these fragments, thus the corresponding partitions 
should not be counted. The number of these invalid 
partitions is obtained by removing one fragment with 
size s max . The other fragments can have any size less 



than or equal to s max and their multiplicity is M — 1. 
The number of invalid partitions is thus Smax N(S — 
M-s max ,M - 1). 

The symmetry relation (Q) can be demonstrated 
graphically using the so-called Ferrers diagram in 
which a size s fragment is representated by a column 
of s dots and a partition by its set of fragments sorted 
in a decreasing order. We complete the Ferrers dia- 
grams with open dots as indicated in figure 0. The 
number of partitions of the integer S into M parts 
with sizes less than or equal to s max is equal to the 
number of ways of arranging the dots in the thin line 
boxe. As the multiplicity is fixed, the bottom row is 
necessarily full. By a 180 deg rotation the open dots 
play the same role as the black dots. The open dots 
partition will be refered to as complementary partition 
of the black dots (which should not be confused with 
the conjugate partition which is obtained by inverting 
M and s max ). Thus, to each multiplicity M partition 
of the integer S corresponds exactly one multiplicity 
M partition of the M (s max + 1) — 5 open dots. 
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FIG. 1. Graphical sketchs of the three partitions 
({3,3,1,1,1}, {3,2,2,1,1}, {2,2,2,2,1}) of the integer 9 into 
5 fragments with size less than or equal to 3 (thin line 
boxes). The bold line boxes contain all the complemen- 
tary partitions ({3,3,3,1,1}, {3,3,2,2,1}, {3,2,2,2,2}) of the 
11 open dots with the same multiplicity and maximum 
size. The hatched dots do not participate in the enumer- 
ation of partitions (due to the multiplicity constraint). 

For example, in the frame of the minimal infor- 
mation model, it can be interesting to know the 
number of partitions containing a given set N : 
(N 3mlD , . . . , Ns) of "large" fragments (i.e. fragments 
with size greater or equal to s m j n ) supplemented by 
"small" fragments : 



2 



N(N) 



(9) 



An alternative method for enumerating the parti- 
tions with constraint on the maximum size consists in 
using the equivalent of equation (|J) which takes the 
following form in this case : 



M 



~ 1 N(S — M, m). 



(10) 



This equation can be applied recursively s max — 1 
times so that the maximum size in the right hand term 
is 1. Using 1 N(S, M) = 1 if S = M and otherwise, 
one obtains : 
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FIG. 2. Ferrers diagrams illustrating the determina- 
tion of the minima and maxima of the sums in equation 
(pq|). The goal is to determine the range of m,2 when 
mo = mi = 4 and S = 11 (thus R2 = 3). The fragments 
are sorted in decreasing order, thus m.2 (s max — 2) > R2, 
that is, the minimum value of m,2 is 2. The number of 
fragments with size greater than k is necessarily lower than 
the number of fragments with size greater than k — 1, thus 
m2 < mi. Furthermore, their are only R2 units left, thus 
m2 < R2- Finally, the maximum value for m,2 is the min- 
imum of mi and R2 (i.e. 3). More generally, the sum for 



rrik runs from Int(- 



+ 1 to Min(m,fc_i, Rk 



Using the same line of thought on the conjugate 
partition, one obtains the following equation : 
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with mo = M and = 5 — J^iLo 1 m «- ^ n ^ n i s 
equation, is the multiplicity of fragments with 
size strictly greater than k. The determination of the 
range for is illustrated in Fig. 0. 



(12) 



with so = 
size of the v 



Si, Si being the 



= Smax and Rk = S 

;ih largest fragment. The same equations 
hold for N(S, M), fixing s max = S. 
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B. Constraint on the minimum size 

The number of partitions of the integer S into M 
fragments with size greater or equal to s m m will be 
noted SmiD N(S, MY In each event, M s m i n units are 
imposed (in Fig. [3| they correspond to the two lower 
rows). The number of partitions only depends on the 

Ms r 
from 1 to M 



S — M s m ; n remaining units, for multiplicities ranging 



M 



to 



lla N(S,M) - J2 N(S-Ms min ,m). (13) 

m— 1 

Following Eq. (§), this expression can be simplified 



,JV(S, Af) = N(S-M(s v 



1),M). (14) 



The same property can be directly deduced by con- 
sideringthe complementary partition (see open dots 
in Fig. ||). 



m 2 = 3 



3 



oo 
oo 



M = 3 



S-Ms n 



M 



una 
oo 
•o 



pool 

ooo 



FIG. 3. Graphical sketchs of all the partitions ({5,2,2}, 
{4,3,2}, {3,3,3}) of the integer 9 into 3 fragments with size 
greater or equal to 2. The two lower rows play no role in 
the counting of the partitions. The bold line box includes 
all the partitions of the integer 3 (i.e. S — M s m i n ) with 
multiplicities less than or equal to 3. 

The total number of partitions is : 

s^N(S) = E S M=7 a N(S-M( Smin -l),M). 



The boundary conditions are : 

N(0,M ^ 1) = and 
JV(0,1) = 1. 



(15) 



(16) 



C. Constraint on the minimum and maximum 
sizes 

When both the minimum and the maximum size 
of the fragments are fixed, the counting of the parti- 
tions is carried out in the same way as previously : 
M s m in units play no role. The number of partitions 
is the same as that of the integer S — M s m i n into 
fragments with size less than or equal to s max — s m i n 
(Fig. ||). Thus, the number of doubly-conditioned 
partitions is obtained as a sum over the number of 
singly-conditioned partitions. 



<N(S,M) = J2 



^N(S-Ms mixi ,m). (17) 



771= 1 



An alternative expression can be obtained by con- 
sidering the complementary partitions (see Fig. ||) : 



-N(S,M)=J2 



^N(Ms max -S,m). (18) 



m — 1 



Applying Eq. (|Tc]) to |t|) and (|T|) one obtains 
respectively : 



™N(S,M) = 
Finally : 



+1 N(S - M(s min + 1), M) , (19) 
- Smi .+i 7V(M(s max + 1) - S, M) . (20) 
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M = 4 

FIG. 4. Graphical sketch of one partition ({5,4,3,2}) of 
the integer 14 into fragments with size included between 
2 and 5. The two lower rows play no role in the partition 
counting. The bold line box contains all the partitions of 
the integer 6 (i.e. S — M s m i n ) into fragments with size 
less than or equal to 3 and multiplicity less than or equal 
to 4. 



III. CONCLUSION 

In this article, we have provided formulae for the 
calculation of the number of partitions with condi- 
tions on the maximum fragment size (Eq. ([|||)). 
with conditions on the minimum fragment size (Eq. 
(|l5|)) and with conditions on both the minimum 
and the maximum fragment size (Eq. ( ^9|j20| )). To 
demonstrate these formulae, the notion of comple- 
mentary partitions was introduced. The constrained 
partition numbers are notably useful in the analy- 
sis of nuclear multifragmentation. Moretto and col- 
laborators Q have introduced an elegant combina- 
torial procedure to isolate rare events correspond- 
ing to the fragmentation of the atomic nucleus in 
a number of nearly equal size IMF (fragments with 
charge greater or equal to Z min ) supplemented by 
light fragments (fragments with charge less than or 
equal to Z m i„ — 1). This procedure requires the eval- 
uation of the number of partitions corresponding to 
a given sum (Z- lm f) of the charges of a given num- 
ber (M) of IMF. This number of partitions is given 
as ZmiiL N(Z imU M) z ^~ 1 N(Z tot - Z imi ). The total 
number of partitions can be evaluated by the follow- 
ing convolution N(S) = J2 S s niln N(s) s ""-~ 1 N (S - s) . 
In a forthcomming article M we will show how the 
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Moretto charge correlation can be calculated explicitly 
in the frame of the minimal information model. More 
generally these formulae are useful in domains where 
the fragment classes (infinite fragments, evaporation 
residues, light particles, intermediate mass fragments, 
liquid and gazeous phases . . . ) are defined with respect 
to their sizes. 
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